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INTERMEDIATE ENERGY PROTON -DEUTERON 
ELASTIC SCATTERING 

By John W. Wilson 
Langley Research Center 

SUMMARY 

> 

A fully symmetrized multiple -scattering series is considered for the description 
of proton-deuteron elastic scattering. An off -shell continuation of the experimentally 
known two-body amplitudes that retains the exchange symmetries required for the calcula- 
tion is presented. The one -boson exchange terms of the two -body amplitudes are evalu- 
ated exactly in this off -shell prescription. The first two terms of the multiple -scattering 
series are calculated explicitly whereas multiple -scattering effects are obtained as mini- 
mum variance estimates from the 146-MeV data of Postma and Wilson. The multiple - 
scattering corrections indeed consist of low -order partial waves as suggested by Sloan 
based on model studies with separable interactions. The Hamada-Johnston wave function 
is shown consistent with the data for internucleon distances greater than about 0. 84 fm. 

INTRODUCTION 

. Although much attention has been given to the theoretical aspects of nucleon- 
deuteron elastic scattering, at intermediate energies, little progress has been made in 
understanding the cross section or polarization at backward angles (refs. 1 to 8). Multiple - 
scattering series have been applied with reasonable success where the first two terms 
(pickup and single -scattering) provide a good qualitative description of the data (refs. 1, 

2, 7, 8, and 9). Quantitatively this success is largely limited to forward scattering which 
is reasonably determined by various approximations of the single -scattering integral 
(refs. 1 and 3 to 9). The usual pickup term is a complicated expression obtained by apply- 
ing exchange symmetry to the direct multiple-scattering series (ref. 10). The "Born" 
approximation is used to simplify this term to a manageable form (ref. 7). It was believed 
for many years that the failure to obtain quantitative agreement near backward scattering 



was due to the inadequacies of the Born approximation to the pickup term (refs. 1, 2, 6, 

7, and 8). At intermediate angles others concluded (based on studies of double -scattering) 
that multiple -scattering corrections are important (refs. 9. and 11 to 13), but the question 
of multiple -scattering effects at backward angles was left unresolved since typical 
approximations prove inadequate at large momentum transfer. The importance of multiple 
scattering was further confused since successive terms of the multiple -scattering series 
are progressively more complicated so that explicit evaluation was made only for the 
simplest of models (for example, using separable interactions) or with many simplifying 
"approximations in calculations with realistic models (refs. 3 and 11 to 15). Other cor- 
Tections were suggested as of possible increasing importance toward intermediate and 

V 

backward angles which include three-body effects (ref. 16), corrections to.the impulse 
approximation (refs. 5, 14, and 17), off-shell effects (refs. 4, 5, 6, 9, and 17), and 
exchange symmetry effects other than pickup (refs. 7, 8,. 15, and 18 to 21). At energies 
above a few hundred Me V, relativistic effects are expected to become: important as well as 
resonance processes. such as those exhibited in the Kerman and Kisslinger wave function 
(ref. 22) and the Craigie and Wilkin double -scattering mechanism (ref. 23). 

A comprehensive treatment of nucleon -deuteron scattering must include an accurate 
description of the scattered nucleon polarization which is a sensitive function of the' phase 
of the amplitude. At intermediate energy, the polarization is observed to have a maximum 
at about 40° followed by a minimum at 100° with a second relative maximum at about 160°, 
all in the center -of -mass system (refs. 18 and 24 to 28). The first maximum and the 
minimum in polarization as a function of scattering angle are reasonably accounted for by 
single scattering and are related to the phases of the nucleon-nucleon amplitudes in a 

t 

sticking -factor approximation (refs. 5 to 8). A second maximum in polarization near 
backward angles appears as interference between the pickup and single -scattering ampli^ 
tudes which is in qualitative agreement with experiment but is far too small in magnitude 
(refs. 7 and 8). 

The backward scattering data are understood in terms of a simple model which is 
valid for intermediate energies. First, it is considered that nucleon exchange which 
arises from applying exchange symmetry to lowest order rearrangement collision (in dis- 
tinction to pickup which is evaluated through Born approximation) is the principal contribu- 
tion near backward scattering. Second, an accurate calculation of the single -scattering 
integral shows it to be almost two orders of magnitude below nucleon exchange (ref. 29) 
and therefore not an important contribution. Finally, an estimate of the magnitude of 
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double -scattering amplitude shows it to be nearly angle independent from 130° to 180° in 
the center of mass and a non -negligible contribution («45 percent). Since the higher-order 
multiple scatterings are expected to be even more isotropic (ref. 14), the higher-order 
multiple scatterings contribute only to an S-wave. Since single scattering becomes 
increasingly small in the backward hemisphere, the model including only an S-wave plus 
nucleon exchange will describe the data. This model has polarization which necessarily 
vanishes at cos e = -0. 8 and -1. 0 with a stationary value between. When the interference 
between the nucleon exchange and S-wave is destructive, a maximum in polarization is 
produced. This model (applicable only in the backward hemisphere) quantitatively explains 
the existing cross-section data and the second maximum in polarization from 77 MeV to 
198 MeV. Although the zero crossing of polarization for the model ( 6 «s 142°) is somewhat 
different from the observed crossing ( 9 «= 135° ± 6°), the addition of a simple S-wave 
appears to be the missing ingredient needed to understand the data. 

The essence of this S-wave is further clarified through the recent model considera- 
tions of Sloan (ref. 30). Sloan considers the three -particle scattering problem for parti- 
cles which interact through a two-body isospin spin -dependent S-wave separable interac- 
tion. The parameters of the interaction are chosen to fit the low-energy nucleon -nucleon 
data. The interesting feature of the model is that it is exactly soluble. Sloan concludes 
on the basis of the model that the multiple -scattering series converges slowly at inter- 
mediate energies but that the slow convergence is restricted only to the S-wave amplitude. 
Furthermore, all partial waves with St > 2 are defined by the nucleon exchange (not a Born 
approximation) plus single -scattering amplitudes with double scattering contributing only 
to the S-wave and P-wave amplitudes. That an S-wave amplitude should result from mul- 
tiple scattering from spherically symmetric two-body interactions is perhaps not too sur- 
prising. The exact details of Sloan's calculations at intermediate energies are questionable 
since the two -body amplitudes in fact show a good deal of angular structure quite aside 
from questions about the correctness of a separable representation. 

It is the purpose of the present paper to examine the adequacy of a multiple - 
scattering description of proton -deuteron scattering at intermediate energy. Although the 
multiple -scattering series is expected to converge slowly, model calculations show that 
the higher-order multiple -scattering terms contribute only to the low-order partial waves 
(ref. 30). The first two terms, neutron exchange and single scattering, are assumed to 
describe the high -order partial waves completely. Relativistic effects and resonance 
processes depend strongly on the depth at which the deuteron is probed and would appear 
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as modification in the higher-order partial wave amplitudes and therefore distinguishable 
from expected multiple -scattering effects. Three -body effects, which contribute to the 
scattering from only a small region of the deuteron interior (ref. 12), would contribute 
to S-wave scattering at intermediate energies, which makes them indistinguishable from 
expected multiple -scattering effects. Estimates of three -body effects show them to be 
small (refs. 12 and 16). In the present work, the deuteron is assumed to couple only to 
the nucleon channel and the internal structure is assumed adequately defined by a non- 
relativistic wave function. The proton -deuteron scattering is represented by a multiple- 
scattering series in which nucleon exchange and single scattering are calculated explicitly. 
The effects of various approximations of the single -scattering integral on the theory is 
considered. The totality of higher -order multiple scatterings is approximated by a phe- 
nomenological S-wave. The extent of the double -scattering contribution will be estimated 
by adding a P-wave to the model. The requirement of high -order phenomenological partial 
waves would be interpreted as evidence of possible resonance processes or relativistic 
effects not included in the present theory. 

SYMBOLS 

a L n ’^L n structure constants of deuteron wave function, n = 1, 2, 3 
D. deuteron four momentum, amu 

_o 

spin-doublet partial-wave amplitude, amu 
E proton laboratory energy, MeV 

-3 

nucleon current matrix element, /3 = S, T, V, A, P, amu 
deuteron mass, amu 
m nucleon mass (=1.0 amu) 

P (8) nucleon polarization, dimensionless 

P^ deuteron D-state, percent 
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p 

p l 

p 2 

Q 

Q £ 


s, t, u 
s 

S,P 



T 


T 

s 

T 

x 

* 0 ! 

y l 


^JLS 
Z 2’ Z 4 


z 


nucleon four momentum, amu 

struck deuteron constituent -four momentum, amu 

spectator deuteron constituent -four momentum, amu 

momentum transfer in center of mass, amu 

_2 

spin -quartet partial -wave amplitude, amu 

2 

Mandelstam variables for two-nucleon system, amu 
Lorentz invariant, amu^ 

-2 

S-wave and P-wave transition amplitudes, amu 
deuteron sticking factor, dimensionless 
transition matrix, amu ^ 

_2 

single -scattering transition matrix, amu 
nucleon exchange transition matrix, amu'^ 

_2 

two-body amplitudes in nucleon spin basis ( a = np, pp), amu 

2x2 representation of spin spherical harmonics with orbital angular 
momentum L, dimensionless 

usual spin spherical harmonic, dimensionless 

uncoupled spin doublet and quartet amplitudes, dimensionless 

off -shell two-body cosine of center -of -mass scattering angle, dimensionless 
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6 

0 


K 


M 


77 


cr(0) 


a 


r a 

0 


2 

nucleon mass difference, off -shell parameter, amu 

proton -deuteron center -of -mass scattering angle, degrees 

deuteron internal momentum, amu 

projection of deuteron spin onto the Z-axis 

deuteron polarization vector, dimensionless 

nucleon -deuteron elastic differential cross section, mb/sr 

Pauli spin vector, dimensionless 

- 3/2 

deuteron momentum space wave function, amu 

Lorentz invariant two-body amplitudes for scattering from protons (a = p) or 
neutrons (a = n) in Current x Current basis, amu^ 


Subscripts: 


E entering state 

F final state 

L deuteron orbital angular -momentum quantum numbers, dimensionless 

£ proton-deuteron orbital angular -momentum quantum numbers, dimensionless 

n neutron 

p proton 
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Three vectors are denoted by above the quantity. Four vectors carry no special 
annotation. Dot products among four vectors b = (b°, b) and c = (c°, c) are as 
b • c = b c - b • c. Dagger (t) denotes Hermitian conjugation. Circumflex (*) denotes 
unit vector. Star (*) denotes complex conjugation. The trace of the matrix A is denoted 
by < A) . 


MULTIPLE SCATTERING AND IMPULSE APPROXIMATION 

The form of the multiple -scattering series is shown in terms of Feynman graphs 
in figure 1 (refs. 12, 19, and 31). Higher-order multiple -scattering terms and three-body 
graphs (ref. IS) are not shown in the figure. The present form of the impulse approxima- 
tion neglects the three -body graphs (ref. 32). 

In the following p and D refer to the four momentum of the proton and deuteron, 
subscripts E and F refer to entering and final states, and tt(m) is the deuteron polariza- 
tion for deuteron spin projection ju along the z-axis in the deuteron rest frame. The 
natural units are chosen in which Planck's constant h and the velocity of light c are unity. 
Energies are measured in atomic mass units, and the difference in proton and neutron 
mass is ignored (m n = mp ~ 1 amu) except in calculating the deuteron binding energy. 

The contributions of the graphs of figure 1 will be used to calculate the transition 
amplitude, and they are obtained in the nonrelativistic limit in a manner outlined by Gross 
(ref. 32). The result is 

T = T X + T g + S + P (!) 

where the four contributions are 
Nucleon exchange: 

T x = ^|B| + K^jcp^TT E , K 

where 


(-* - 


<P hr 


F’ k F 


( 2 ) 


_ - l n 
C E,F " P E,F " 2 E,F 


B = M D - m p - m n 
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Single scattering: 


where 


T S=J I I ^(s,t,u, 

a=p,n /3=S,T,V,A,P 




1< e = k + Q/4 
K-p = K - Q/4 

Pl E,F = i° E ’ F+ * E > F 

Higher order: 


S - D 0 Z 2 + Qo z 4 (^F’ ^e) (4a) 

P=Dl cos 0 2^2 ) + Ql cos ® z 4 (^F’ ^e) (4b) 

with the following explanation of the functions and variables. 

The first graph on the left in figure 1 is the well-known nucleon exchange term 
which results as a direct application of exchange symmetry in the lowest -order rearrange- 
ment graph (order refers to the number of proper vertices). This term contributes to the 
scattering in the elastic channel only because the deuterons arising from direct and rear- 
rangement collisions are experimentally indistinguishable. This amplitude has been pre- 
viously derived through Feynman analysis (ref. 33), and is usually expressed in terms of 
the Blankenbeckler-Cook amplitude (refs. 32 to 35), which has in turn been related to the 
nonrelativistic wave function by Gross (ref. 32). Thus, the result is equation (2) in the 
nonrelativistic limit with <^(tT, /T) representing the wave function (to be discussed more 
fully) of a deuteron with polarization if and relative internal momentum if. The pickup 
term of potential theory arises from applying exchange symmetry to the direct multiple- 
scattering series terms (ref. 10). A term equal to the nucleon exchange term is obtained 
as a Born approximation to the pickup amplitude (refs. 7, 8, and 14). 

The next two graphs in figure 1 represent the scattering of the incident projectile by 
a single constituent of deuterium. All vertices are properly symmetrized under exchange 
of identical nucleons. The single -scattering amplitude has been previously derived 
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through Feynman analysis (refs. 11, 12, and 19). The nonrelativisitic reduction follows 
closely the work of Gross (ref. 32). The nucleon -nucleon scattering amplitude is retained 
in fully relativistic form. The two-body amplitudes are represented as (refs. 36 and 37): 


a(Pl F -P2 F ;Pl E -P2 E )= I T s( s ' t ’ u > 6 >jp(Pl F ,Pl E ) X i' i (P2 F .P2 E ) 


(3=S,T,V,A,P 


where p is summed over the five /3-decay type currents; Scalar, Tensor, Vector, Axial- 
vector, and Pseudo -scalar. A similar Current x Current representation appears in ref- 
erence 19 by Chahoud and others. The s, t, and u variables are the usual Mandelstam 
variables and 5 is the difference of mass squared of the struck constituent in the entering 
and final state. The two -body amplitudes are the symmetrized nucleon-nucleon ampli- 
tudes. This symmetry is expressed through 

T*(s,t,u,6) = (-l) 1 ^ 

a 

where I = 0, 1 is the isospin state and is the Fierz transformation. Clearly any off- 
shell approximation should preserve this symmetry (that is, the Pauli principle). The 
continuation of experimental data to off -shell values is given in a later section (see also 
ref. 29). A term similar to single scattering results from applying exchange symmetry 
to rearrangement processes in potential theory (refs. 30 and 38). 

The totality of higher-order terms is represented by the S-wave of equation (4a). 
This representation is in accord with the work by Sloan (ref. 30) who showed through 
model calculations that triple and higher-order terms in the multiple -scattering series 
converge slowly and contribute to an S-wave. The double scattering is suggested by Sloan 
to contribute most to the S-wave but with a few -percent contribution to the P-wave of 
equation (4b). The P-wave contributions coupling the doublet and quartet spin states are 
neglected. The two uncoupled spin transition amplitudes are given by 


Doublet 
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where o is the ordinary Pauli spin operator. The and of equations (4a) and (4b) are 
complex functions of center -of -mass energy which are presently uncalculated. They will 
be determined so as to give best fits to the experimental data. 

Calculation of proton -deuteron elastic scattering in an impulse approximation 
requires, as inputs, information concerning the two-nucleon system, that is, the deuteron 
wave function and two-body scattering amplitudes (compare eqs. (2) and (3)). For the 
energy range considered here, these inputs are relatively well known from more than 
three decades of experimental work. Much of this experimental information is summa- 
rized in recent phenomenological potentials (refs. 39 to 42) and phase-shift analysis 
(refs. 43 to 45). 


DEUTERON WAVE FUNCTION 

The most general form of the deuteron wave function is implied from the following 
invariant operations: translation, rotation, space inversion, time inversion, and nucleon 
exchange (Pauli). From translational invariance, it is concluded that the wave function of 
a deuteron of momentum D, polarization 7, projected onto a two- nucleon state has the 
representation 

%i' 2 (* 1 ’* 2 ’ 5,5r ) = ( 27r )" 3/2e ’ 1X ' 5 

with 

Va^’* 1 ’* 2 ) = (27r) " 3/2 I d3/<e ’ ir %^ 2 ( 7 ’Pi’P2) 

and 

x = |(x 1 + x 2 ) 

7 = x x - x 2 

« = |(Pl'P2) 

where x-^ and x 2 are the position vectors of the two nucleons and and i^ 2 are their spin 
projections. 


(7) 

( 8 ) 
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Considered as a 2 x 2 matrix, 


, *' V2 ( 7 ’Pl’P2)J = ‘K’.Pl.^K 1 


( 9 ) 


with 


c = 


0 -1 

1 0 


It is found that cp(ir, Pj, p 2 ) has the following form in the deuteron rest frame (pj = -p 2 = k) 
due to the remaining invariances: 


<p{*,k) = <p 0 {k)y 0 (1?) - (p 2 (l<) y 2 (7,k) 


( 10 ) 


with 


&Y 0 (f) = n • ct J\/2 

i/4^Y2(tT, k) = 1 5T • k k • a - ^ It • o 


The matrices Yq and Y 2 are related to the ordinary spin spherical harmonics by 


l -{^H V"2 -*£/*) 

v l v 2 1 z 

I {y s [w,i}-‘) sS = i. (!: 
v \ v 2 1 z 


(11) 


(12) 


where Xj/2 = a ’ X-l/2 = P is in agreement with the representations of Gross (ref. 32) 
and Gourdin and others (ref. 34). It is concluded from equations (10), (11), and (12) that 

<P 0 (*> = /i io j o (Kr) u(r) r dr 

<P 2 (x) = ‘/I j 2 (/cr) w(r) r * 

where u(r) and w(r) are the ordinary S and D state radial wave functions. 
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Various solutions for u(r) and w(r) are offered for different phenomenological 
potentials (refs. 40, 41, 42, and 46). The solutions to the Hamada- Johnston potential 
found with variational methods by Humberston and Wallace (ref. 46), which are presented 
in a convenient analytic form, will be used here. The Fourier transform of their analytic 
wave functions (ref. 46) can be written as sums of exponential integrals of complex argu- 
ment as used in the present calculations. 

TWO-BODY TRANSITION AMPLITUDES 

The second input to the impulse approximation is the nucleon -nucleon scattering 
amplitude. Although this amplitude is not known from fundamentals, it has been the object 
of much experimental work and is reasonably well known to a few hundred MeV with 
uncertainties entering at higher energies due to lack of experimental data and knowledge 
of inelastic processes (ref. 44). 

The on-shell two-body amplitudes r^(s,t) are calculated from the Livermore 
phase shifts (ref. 45) but must be continued off -shell before the integral, equation (3), can 
be evaluated. 

The invariant off-shell amplitudes are taken as functions of the three Mandelstam 
variables s, t, and u with a fourth variable 6 as the mass difference of the struck nucleon 
in the entering and final state 

S = ( P E + %) 2 
t = ( p F ' Pe/ 

U = ( p i F " p e ) 2 

. _ ... 6 =V' P lE 2 

In general, the values of the four invariants are independent and 
s + t + u ^ 4m 2 
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The on -shell -like variables s and S 0 defined by 
s + t + u s 4m^ 

6o = ° 

are introduced, and the invariant amplitude is expanded in a Taylor series in the neighbor 
hood of the on -shell region 

Tp(s, t, u, 6) = t“(s, t, u, 6q) + (s - s)j- t“(s, t, u, 6 0 ) + (6 - 6 0 )^r t“(s, t, u, 6 0 ) + . . . 

For the present, only the first term of this series is retained which is the experimentally 
known amplitude. This off -shell continuation evaluates exactly the one -boson exchange 
parts of the amplitudes and preserves the Pauli principle in the continuation (for more 
detail see ref. 29). The errors in neglecting the second term were shown to be small 
(at least for quasi-elastic scattering) at intermediate energy through explicit evaluation 
(ref. 17) while the third term must vanish if time reversal invariance is to be preserved. 

EVALUATION OF THE SINGLE -SCATTERING INTEGRAL 

The principal assumption normally used in developing approximations to single 
scattering, equation (3), is that the product of two -body amplitudes and spin spherical 
harmonics is a slowly varying function of internal momentum compared to the radial - 
wave functions. Taylor series expansions of two -body amplitudes and spin spherical 
harmonics are then made with only the first term retained, leaving an integral over radial 
wave functions (the sticking factors). Various approximations are obtained by the choice 
of expansion point. The only approximation of this type considered here is the Kottler- 
Kowalski prescription (ref. 7). An alternative is to examine the radial-wave function 
products which have many relative extrema. The two-body amplitudes and spin spherical 
harmonics are then expanded about each extremum and multiplied by appropriate weight 
determined from the deuteron wave functions (refs. 29 and 47). This procedure results 
in an approximation which combines sums of two-body amplitudes and spin spherical 
harmonics evaluated at points of principle contribution to the integral, equation (3). A 
third method is to perform a numerical integration in which a convergence test is used as 
a standard for adequacy of approximation (refs. 29 and 47). 
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The implementation of the above approximations of the integral in equation (3) is 
aided by taking the deuteron wave function as a sum of Gaussians 


“ Z a L n exp (■ b L„' f2 ) 


The coefficients given in table I are fits to the Hamada -Johnston wave functions (ref. 46) 
and are reasonable fits to 0. 28 fm. It is seen from figure 2 that the S -state is least accu- 
rately approximated. As a check on the effect of these errors on the integral in equa- 
tion (3), the sticking factors for the S-state are compared in table n to the analytic 
sticking factors derived from the Humberston-Wallace L = 4, N = 2 analytic wave func- 
tions (ref. 46) and the values computed by Glendenning and Kramer for their potential 
number 9 (ref. 48). 

The three approximations of the integral in equation (3) are represented as follows: 
Let the product of two -body amplitudes and spin spherical harmonics be combined into the 
symbol 

T (s, t, u, 6 )l l f = ^ X T ^ S ’ U ’ ^ ^( P F> p e) 

a=p,n 0=S,T,V,A,P 


X <Yl E (* e > k e ) ^Lf^F’ k f) j ( P 1 F ,P 1 E | 

so that the single -scattering integral in equation (3) becomes 

Ts= E I T(s ’ t)U ’ 6 )L E L F < ^L E ( K E) <p L F ( K F) d3,< 


l e l f 


■ X 


l e l f 


a T a T 
L E n L Fm 


exp^-Q!(l - p 2 )Q 2 /l6 J r(s,t,u, 6 ) LeLf exp^a^K - 


where the last equality follows from use of the Gaussian wave functions and 


«= b L F + b L E 
* m ^n 


P = [ b L^ ' b L 


F m E n 
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The following approximations may now be written: 

Kottler -Kowalski (ref. 7) 

T s ~ ^ T ( s ,t, u , <5) Le L f s(Q)l e L f ( 18 ) 

l E l F 

where the sticking factor is given by 


l E l F 




- 2 a L E a L F (!)<*p[“(l - /5 2 )Q 2 / 16 ] 
n, m n m 


(19) 


Asymptotic method 


T s* Z a L E a L F (5)^|-«( 1 -^ 2 )Q 2 / l6T < s » t » u * 6 )L B L F 
n r m L ' 


Le l F 

n m 




( 20 ) 


Numerical integration 

The variables chosen for the numerical integration are for a cylindrical system 
with z along Q. This choice reduces equation (3) to a form for direct application of Gauss 
integration by using Hermite, Laguerre, and Legendre polynomials for components of /T 
along Q, for components perpendicular to Q, and for the angular variable, respectively. 

The magnitude squared of the single -scattering integral for the three approxima- 
tions is shown in figure 3. There is a small error for the Kottler -Kowalski and 
Asymptotic -method approximations in the forward hemisphere due to inadequate approxi- 
mation of the integral over the D-state spin spherical harmonic. The principal errors in 
the backward hemisphere are due to the rapid variation of the two-body amplitudes which 
arise from exchange symmetry terms in the two-body scattering (refs. 29 and 47). The 
main result from accurate evaluation of single scattering is the continued decrease of 
the amplitude in the backward hemisphere. The errors of low-order approximations 
have been incorrectly interpreted by others as the mechanism producing the backward 
scattering peak at higher energies (refs. 19 and 20). 
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ANALYSIS OF EXPERIMENTAL DATA 


The proton -deuteron elastic scattering amplitude (eq. (1)) is related to the dif- 
ferential cross section o(6) and the proton polarization produced by the scattering P(0) 
by 


O(e) = (2tt) 4 C 



( 21 ) 


P(0) = ^ (rt a z T^> / ^ <T* t) (22) 

Me^f J ^e^F 

where <( ^ denotes trace over nucleon spin indices and C = 10(0. 2101)^ is a conversion 
-2 

factor from amu to mb. The theory includes parameters (principally the S and P wave 
contributions of higher-order multiple scattering) which will be adjusted to obtain best 
agreement with the data. These parameters are estimated by minimizing 



where P is the parameter estimate, the quoted experimental uncertainties (Acr eX p and 
AP eX p) are chosen as weights, and the uncertainty due to finite detector size is ignored. 
The factors N a and Np allow for renormalization of experimental data and relate to the 
uncertainty of initial beam intensities (unpolarized and polarized) used in the experiment 
which are often known to no better than 10 percent. 


Assuming the experimental errors (for fixed normalizations) to be normally dis- 
tributed with standard deviation given by the quoted experimental uncertainty makes the 
minimum F(P) to be bounded by a chi-square variable (actually, F(P) ~ x^) with the num- 
ber of degrees of freedom equal to the number of values of cross section and polarization 
combined. For N x (number of data points) large, the test for goodness of fit requires 


16 



N x - 2/2% § F(P) g n x + 2 / 2 


for the fit to be acceptable. Standard minimization procedures were used to find the 
minimum. 


RESULTS AND DISCUSSION 

The principal parameters in the present model are the S-wave and P-wave ampli- 
tudes of equation (1). In addition, the D-state percentage is taken as a parameter in the 
model in an attempt to extract it from the experimental data. The shape of the radial 
wave function is taken to be that of Hamada and Johnston (ref. 46). The parameters are 
estimated from the 146-MeV cross-section and polarization data of Postma and Wilson 
(ref. 18). The polarization renormalization is taken as Np = 0.93 in accord with refer- 
ences 25 and 49. The cross-section renormalization N a was taken as a parameter and 
is to be compared to the latest quoted value N a = 0.92 (ref. 26). 

The data analysis was first guided by the suggestion of Sloan (ref. 30) that higher- 
order multiple scattering is predominantly S-wave with double scattering contributing a 
small P-wave (few percent) amplitude. The initial fits ignored the P-wave completely. 

The data were fit using only the S-wave of equations (1) and (21) to (23) for the three 
approximations of the single -scattering, integral. The resulting fitted parameters are 
given in table III with the resulting curves in figures 4 through 6. The experimental 
values (rectangles in figures) plotted in the figures are not those of the originators 
(ref. 18) but have been renormalized by N ff = 0.92 (ref. 26) and Np = 0.93 (refs. 25 and 49). 
It is seen from the values of F(P) obtained that the more accurate approximations of single 
scattering provide progressively worse agreement with the data. This fact can also be 
observed in the visual comparisons, especially at backward angles, in figures 4 through 6. 
The deterioration of the fits for the more accurate approximations occurs since backward 
single scattering interferes destructively with the nucleon exchange amplitude. A more 
accurate approximation of single scattering diminishes its magnitude (see fig. 3), hence 
reducing the interference term and thus providing worse fits to the data. The values of 
F(P) for each of these three fits are outside the acceptable range given by 60 ± 22. The 
assumption that higher-order multiple scattering is mostly S-wave must be rejected. 

The fact that higher-order multiple scatterings are not predominantly S-wave as 
predicted by Sloan (ref. 30) is not too surprising since his predictions are based on the 
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assumption that the two -body amplitudes are isotropic and separable. Such a representa- 
tion of the two -body amplitudes is not realized in nature. It is believed that the P-wave 
contribution from double scattering predicted by Sloan is underestimated and reasonable 
fits should be obtained by adding a non-negligible P-wave to the model. This possibility is 
contingent on relativistic effects and resonance processes being negligible at this energy 
since their first-order effects are expected to modify high -order partial-wave amplitudes. 

The curves in figure 7 show the effect of adding the uncoupled spin P-wave to the 
formalism. It is seen in table III that the minimum of the fitting function drops to a value 
which indicates that the average data point is about 1 o distance from the theoretical curve. 
The parameters were fit again using the numerical integration for the single -scattering 
integral with P-wave included but with the D-state percentage restricted to 6.93 ± 1. 0 and 
6. 93 ± 1. 5, and the respective results are shown in the last two lines of table III. The 
extreme (Pp> = 5. 93) of these last two results is compared in figure 8 with the experimental 
data. Comparison of the results in figures 7(a) and 8(a) shows the major adjustment to be 
in the cross sections between 140° and 170° scattering angles. A modest improvement in 
cross section between 40° and 80° for the higher Pq is also noted. It is impossible to 
distinguish on the basis of a chi-square test between these last three fits. 

Examination of table III and the visual comparison of figures 6 through 8 show that 
the addition of a large P-wave contribution to the model is unavoidable. A rather large 
range of D-state percentage (3.93 = Pj) ^ 5.93) is shown consistent with the data. That 
the D-state percentage is poorly determined is due to correlations among the fitted param- 
eters and the fact that the D-state contribution is a second-order effect at this energy. 

More accurate cross sections in the range of 40° to 80° and 140° to 170° could possibly 
help to reduce the range of Pq over which good fits can be obtained (compare figs. 7(a) 
and 8(a)). An alternative may be to calculate explicitly the main contribution of the S- and 
especially the P-wave amplitudes (that is, double scattering) and thereby eliminate the 
principle part of the correlations in the fits. However, such a calculation may not be 
possible with full spin dependencies on present-day computers unless some innovative 
approximations are found. 

Finally, the extent of the multiple -scattering effects are summarized in figure 9 
by comparing the partial cross sections of nucleon exchange plus single scattering with 
the multiple -scattering contribution estimated from the data (that is, S + P). It is seen 
from figure 9 that the multiple scattering has two broad maxima at forward and backward 
angles. The cross section is dominated by multiple scattering in the region 110° to 140°. 
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Presumably, the angular structure of the multiple -scattering amplitude would be given by 
double scattering. It would be interesting to see whether a meaningful calculation of 
double scattering would give such results. 

CONCLUDING REMARKS 

The present belief in the correctness of a multiple -scattering series formulation 
of proton -deuteron scattering lies in the fact that the first few terms plus low -order par- 
tial waves provide a model in good agreement with the data. The fact that high -order 
phenomenological partial waves are not required indicates that relativistic corrections 
and resonance process are not important at the present energy. The low -order partial 
waves are understood to be higher -order multiple -scattering terms on the basis of model 
calculations of Queen and especially Sloan. The large phenomenological P-wave contribu- 
tion to the amplitude in comparison to the few percent predicted by Sloan is believed due 
to the much stronger angular variation of the real two-body amplitudes as opposed to the 
S-wave two-body interaction in Sloan's model. The clarity of current understanding lies 
in the present structure of the multiple -scattering formalism in which nucleon exchange, 
fully symmetrized single scattering, etc. , arise in a natural way from Feynman analysis 
of the three -body amplitude. In distinction, the pickup term is usually obtained as a Born 
approximation to exchange symmetry terms of the usual direct multiple -scattering ampli- 
tude. This procedure led to the widely accepted conclusion that inadequacies of the theory 
were indeed a problem of the Born approximation. The present result eliminates this 
possibility and clearly indicates higher-order multiple scattering as the necessary ingre- 
dient to complete the theory provided, of course, that terms are properly symmetrized. 

It has been shown that the D-state is poorly determined from the presently exam- 
ined experimental data. This poor determination occurs because the range of momentum 
over which the deuteron is presently probed is mostly the region where the S- state con- 
tributions are most important (Internal momentum = 0.25 amu). Those regions where the 
D-state is most effective in nucleon exchange and single scattering are also the regions 
where multiple -scattering effects are very important. Although there is modest improve- 
ment in cross section in the angular range 40° to 80° in going to higher values of D-state 
percentage, both values of D-state percentage are seen to be within the experimental- 
error limits. A positive conclusion from the fits is that the Hamada- Johnston S -state 
wave function is shown consistent with the data to distances beyond about 0. 84 fm. 


19 



Although it is gratifying to see that a multiple -scattering theory of proton -deuteron 
scattering is nearly at hand, at least at intermediate energy, it still remains to be demon- 
strated that inclusion of explicitly calculated higher -order terms will bring about improve- 
ment of the theory. Even if improved calculations are made, it currently appears that 
little additional information about the deuteron wave function will be learned at the pres- 
ently considered energy. Perhaps the most important result from the present study is 
an estimate based on the Postma and Wilson experiments of the multiple -scattering effects 
to which future calculations will hopefully converge. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 24, 1973. 
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TABLE I.- FITTED STRUCTURE PARAMETERS FOR GAUSSIAN APPROXIMATION 

OF HAMADA-JOHNSTON WAVE FUNCTION 


Normalization 



dq = 


i; 


= 2 a L n eXP (' b L n q2 ) 

n 


n 

-3/2 

, amu 
n 

k -2 

Dq , amu 
n 

-3/2 

ag , amu 
n 

k -2 

t >2 , amu 

n 

i 

92.41770 

204. 15580 

-1. 683110 

11. 27690 

2 

12.74490 

25.77230 

-. 971535 

2. 85898 

3 

-.35298 

1. 29835 

2. 659846 

465. 61200 


TABLE U. - COMPARISON OF STICKING FACTORS OF APPROXIMATE GAUSSIAN 
WAVE FUNCTIONS WITH ANALYTIC HAMADA-JOHNSTON SOLUTION 
(REF. 46) AND THAT OF GLENDENNING AND KRAMER (REF. 48) 


|Q|, fm” 1 

S-state sticking factor, S(Q)q q 

H-W* 

Approximate 

Gaussian 

G-K f 

0.0 

0. 930 

0.928 

0.926 

.5 

.797 

.790 

.795 

1.0 

.555 

.526 

.557 

1. 5 

.356 

.323 

.361 

2. 0 

. 219 

. 201 

. 226 

3.0 

.071 

.064 

.078 

4.0 

.011 

.007 

.016 

5.0 

-. 010 

-.010 

c— 

o 

o 

1 

6.0 

-. 016 

-.012 

-.014 


*Humberston- Wallace L = 4, N = 2 analytic solution. 
^Glendenning and Kramer values for potential number 9. 
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Figure 1.- Fully symmetrized multiple -scattering series 
in terms of renormalized fermion lines and proper 
vertex parts. 



k, amu 


Figure 2.- Comparison of approximate Gaussian wave function with 
the Hamada- Johnston wave function obtained by Humberston and 
Wallace (ref. 46). 





Figure 3.- Modulus squared of single-scattering integral as function 
of center-of-mass scattering angle for 146-MeV proton laboratory 
energy. 
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Figure 6.- Comparison of S-wave fits for numerical integration of single scattering. 
Best estimate of D-state percentage is 5.16. 
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Figure 9.- Comparison of partial cross section of nucleon exchange plus single 
scattering to that produced by multiple scattering as estimated from experi- 
mental data. 
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